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ABSTRACT 

In  thu  wort,  we  diecuee  the  eottttioa  of  the  differential  equauone  in  the  cylindrical  coordinate  system 
for  the  cue  of  e  coupled  CE  rcacrioa  medienmm  at  a  finite  duh  electrode.  We  use  the  properuee  of 
disoonnnuoue  intefralt  of  Bernal  fuactiana  to  treat  the  mixed  boundary  condition!  at  the  electrode 
surface  and  the  surroundint  insulator.  The  reeulta  are  compared  to  previous  analyses  and  we  also  make 
comparisons  of  various  estimates  of  the  chronoampcromemc  response  (or  constant  surface  concentreuon 
conditions.  ■ 


INTRODUCTION 

The  application  of  microdisk  electrodes  to  the  investigation  of  the  kinetics  i 
chemical  reactions  in  solution  coupled  to  the  electrode  reactions  (CE  ECE  r.d 
D1SP  reactions)  has  been  investigated  recently  [1,2].  It  was  shown  that  the  dimer, 
sion s  of  the  electrode  become  a  parameter  of  the  investigation,  these  dimensions  r. 
effect  probing  the  reaction  layer  at  the  electrode  surface.  An  interesting  aspeo: 
these  measurements  is  that  it  becomes  possible  to  differentiate  between  alterru; . 
reaction  paths  such  as  the  ECE  and  DISP1  mechanism  by  combining  measuremer.- 
using  microelectrodes  with  those  using  conventional  planar  electrodes;  rate  para:-  r 
ters  can  be  derived  from  “working  curves”,  e.g.  of  the  effective  number  of  elec;: 
transferred  as  a  function  of  the  inverse  of  the  radius  of  the  electrode. 
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\n  approximate  method  of  data  analysis  was  used  in  these  investigations.  For 
example,  for  the  CE  mechanism 

A«B  (*> 

ks 

B-i^C  (u) 


the  flux  of  B  to  the  surface  was  obtained  from  the  solution  of  the  differential 
equauon  in  the  spherical  coordinate  system  in  the  steady  state 


2D  dc 
r  dr 


+■  k\  -  k^c  “  0 


0) 


rather  than  that  for  the  cylindrical  coordinate  system  appropriate  to  the  microdisk 
electrode 


3c  _  3Jc 
—  —  D— 
3  i  3  r1 


r  Or  dzi  1  * 


(2) 


in  the  steady  state.  The  solutions  obtained  for  a  microsphere  were  then  corrected  to 
those  for  a  microdisk  by  assuming  that  the  same  scale  factor  holds  in  these 
experiments  as  for  the  measurements  of  the  voltammetnc  curve  for  a  reversible 
electrode  reaction,  viz. 


^diafc  * 


(3) 


where  and  are  the  relevant  electrode  radii 

We  have  shown  recently  that  it  is  possible  to  develop  systematic  algebraic 
analyses  of  a  variety  of  the  usual  types  of  electrochemical  experiments  (chronopo- 
tentiometry  [3],  chronoamperometry  [4],  ac  impedance  [3D  by  appropriate  solutions 
of  the  differential  equations  in  the  cylindrical  coordinate  system.  In  the  derivation 
of  these  solutions  we  have  made  use  of  the  properties  of  discontinuous  integrals  of 
Bessel  functions  to  take  account  of  the  mixed  boundary  conditions  at  the  electrode 
surface,  viz.  prescribed  concentration  or  flux  over  the  electrode  surface,  zero  flux 
over  the  surrounding  insulator  surface.  In  this  paper  we  explore  the  application  of 
the  same  approach  to  the  analysis  of  the  kinetics  of  chemical  reactions  coupled  to 
electrode  reactions  using  the  CE  process  (i)  and  (ii)  to  illustrate  the  method  of 
calculation. 


THEORETICAL  CONSIDERATIONS 


In  the  simplest  method  of  investigating  the  kinetics  of  the  coupled  chemical 
reactions  (i)  and  (ii),  a  sufficiently  high  overpotential  is  applied  to  the  electrode 
surface  to  reduce  the  surface  concentration  of  B  to  zero.  The  rate  of  the  reaction 
then  becomes  controlled  by  the  kinetics  of  reaction  (i);  species  A  is  assumed  to  be 
present  at  sufficiently  high  concentrations  that  we  need  only  consider  the  concentra- 


lion  distribution  of  species  B.  eqn.  (2).  It  should  be  noted  that  it  is  relatively 
straightforward  to  achieve  such  conditions  with  microelectrodes  (and  of  pseudo- 
first-order  reacuon  conditions  in  second-order  reactions)  in  view  of  the  high  rates  of 
steady  state  mass  transfer  to  the  electrode  surface  (1,2]. 

We  therefore  require  the  solution  of  eqn.  (2)  subject  to  the  condition 

c  — /•—  oo.  a  —  ao.  all  t  (4) 


and  with 


c  —  0  0<r<a.  r  —  Q,  t>0  (51 

3c/3r«0  /■>  a,  r-O,  r>0  (6) 

We  have  pointed  out  that  the  application  of  condition  (S)  to  the  solution  of  the  mass 
transfer  problem  for  the  chronoamperomethc  response  is  difficult  and  these  difficul¬ 
ties  also  apply  to  the  present  case.  Instead  we  consider  that  the  electrode  is 
subjected  to  a  constant  uniform  flux  -  Q,  Le, 

2>3c/3?--g  0<r<a,  z-0,  r>0  (7) 


and  we  evaluate  the  average  concentration  of  B  over  the  electrode  surface  and  then 
determine  the  strength  of  Q  required  to  make  this  average  concentration  zero.  It 
should  be  noted  that  the  constant  surface  flux  condition  is  in  fact  likely  to  be  a  very 
good  approximation  for  the  investigation  of  coupled  chemical  reactions  since  the 
rate  will  be  controlled  by  step  (i).  Provided  the  concentration  of  A  is  uniform  over 
the  surface,  Q  will  be  uniform  over  the  surface. 

Laplace  transformation  of  eqn.  (2)  subject  to  condition  (4)  gives 


3 1  3c 

3,3  +  7  3 r 


(8) 


where  q  -  (s/D)l/1. 

The  solution  of  eqn.  (8)  consists  of  the  Particular  Integral  kt/k2s  and  the  Comple¬ 
mentary  Function  which  must  be  found  from 


d^cF  +  2.  d^cr 
3,3  r  dr 


7>)?CF“° 


(9) 


By  analogy  to  the  discussion  of  the  previous  examples  of  electrochemical  experi¬ 
ments  at  microdisk  electrodes  [3-5],  we  note  that  the  substitution 


cCF-u 


converts  eqn.  (9)  into 


(10) 
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with 


/  K* 


-(*=♦£) 

^  A*  1 

(12) 

^)«p[/(^.?.§)rj2y0(ar)da 

(13) 

The  solution  is  then 


*  kyD)  be‘*“- »  «-y  *•  ™„o™ 


dc 

dr 


-0. 


/■>  a,  z-0 


«-  -g(Q  n 

dr  D —  0<r<a,z-0 


We  again  make  use  of  the  discontinuous  integral 
-00 

/  Jo(ar)Mat>)  da- 

'0 


l/2a 

ll/a 

and  with  the  particular  form 

[/(*.  *,£)]' -a 

We  obtain 


r>  a 
r  —  a 
r<a 


■m  *V‘2^reXP['(aJ  +  ?,  +  ^)',/J]^)A(aa) 

v-  da 


<-«-r 

and,  at  z  -  0, 

A.  g(*)q  /•”. ,  , 

kji  £>  ;0  /o(ar)y,(aa)- 


da 


We  derive  the  average  surface  concentration  over  the  surface  of  the  disk 

'd'“« 


A  _  2Q(*)  f°°r  .  / 


(“2  +  ^  +  ^ 


Z> 


da 


a(ai+*i+£r 


(14) 

(15) 

(16) 

(H) 


(18) 


(19) 


.(20) 
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and.  on  setting  the  average  concentration  at  the  surface  equal  to  zero,  we  obtain 
k,D/2kzs 


Q(s) 


da 

/  [  /,(  aa  ))  - 

0  a(  a:  +  qz  +  k^/D) 


If  we  define 

D'  'a 


(s*k.y 

eqn.  (21)  can  be  written 


£>(*)■ 


hl.ai 

2k,  “1 

f  r  +  lcj) 

D  J 

1/2 

1 

»i:~ 

8 

[H 

is+-ki) 

n 

?  dfi 

{  D  j 

(21) 


(22) 


(23) 


We  therefore  obtain  a  family  of  curves  in  the  r-domain  of  Q(s)  as  a  function  of 
(a2s/D)  each  determined  by  the  given  value  of  the  parameter  (aJifc 2/ D).  Numerical 
inversion  gives  the  responses  in  the  r -domain.  Here  we  restrict  attention  to  the 
steady  state  limit 


(24) 


It  is  readily  confirmed  that  this  result  also  follows  directly  from  eqn.  (2)  in  the 
steady  state  provided  we  use  the  substitutions 

(25) 


c  «*  v  exp 


and 

$  -  Dx/ia/k\n 


(26) 


Equation  (24)  is  identical  to  the  chronoamperometric  response  [4]  provided  we 
replace  (c*  -  cAv)  by  and  sl/l,  by  kV1,  i.e.  we  obtain  - 


(2") 


The  result  of  the  approximate  method  of  analysis  (see  above)  can  be  written  in  the 
form 


(2Si 


* 

~> 

L 
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TABLE  1 


Comparisons  of  the  various  results  for  the  CE  flux  and  diffusion  limited  current  flux 


a:k , 

Exact  eqn. 
C7) 

Equation 

1 28) 

Equation 

(29) 

(eqn.  29) -(eqn.  27) 

(SS)‘  -  ( eqn. 27) 

D 

(eqn.  27) 

(eqn.  27) 

0.01 

23701 

25664 

23747 

0.001943 

0.000834 

0.015 

10564 

11451 

10595 

0002901 

0.001248 

0.02 

j960.1 

6466.2 

5983.0 

0.003849 

0.001660 

0.03 

2664.1 

2896.1 

2679.1 

0.005717 

0002476 

0.04- 

1507.5 

1641.5 

1518.9 

0.007547 

0.003284 

0.05 

970.42 

1058.6 

979  49 

0.009341 

0.004082 

o.oc 

677.82 

740.69 

685.34 

0.011098 

0.004872 

0.07 

500.87 

548.26 

50709 

0.012821 

0.005653 

o.os 

385.69 

422.89 

39109 

0.014508 

0.006425 

0.09 

306.50 

336.60 

311.43 

0.016162 

0.007189 

0.1 

249.69 

274.65 

254.12 

0.017781 

0.007945 

0J 

66.028 

73.662 

68.151 

0.032251 

0.015056 

0.3 

30.987 

34.961 

32.341 

0.043929 

0.021421 

0.4 

18.374 

20.915 

19.333 

0.053256 

0.027117 

0.5 

11376 

14.186 

13.126 

0.060606 

0032211 

0.6 

90306 

10.407 

9.6292 

0.066293 

0.036765 

0.7 

6.9609 

8.0540 

7.4322 

0.070585 

0.040133 

0.8 

5.5832 

6.4789 

5.9947 

0.073709 

0.044462 

0.9 

4.6150 

5.3660 

49650 

0.071854 

0.047696 

1.0 

3.9053 

4.5463 

40067 

0.077185 

0.050573 

2.0 

14186 

1.6366 

13143 

0.067456 

0.065799 

3.0 

0.83920 

0.94961 

0.87862 

0.046943 

0.067957 

40 

0.59278 

0.65915 

0.60990 

0.028887 

0.065204 

5.0 

0.43767 

0.50186 

0.46436 

0.014619 

0.060779 

6.0 

0.37258 

0.40407 

037387 

0.003464 

0.056015 

7.0 

0.31414 

033761 

031245 

-0.005360 

0.051451 

8.0 

0.27152 

0J8979 

006813 

-0.012472 

0.047287 

90 

003908 

005366 

0.234/0 

-0.018303 

0.043569 

10.0 

0.21356 

002546 

0.20862 

-0.023153 

0.040282 

12.0 

0.17599 

0.18435 

0.17058 

-0.030741 

0.034831 

*  Simulation  result  from  Shoup  and  Sabo,  reference  (7). 


The  result,  however,  applies  lor  constant  surface  concern! alien  -hcicc;  eqn.  (fP) 
assumes  a  constant  surface  flux.  The  correction  to  the  latter  condition  will  be  of  the 
order  of  the  ratio  of  the  mass  transfer  coefficients  for  constant  flux  and  constant 
concentration,  3rr2/3Z  i.e.  we  predict  the  relation 


„  M  3 ir/  D  W,  7t  I  a2k1\'/i\ 

Q~  2  '  <  (a2*J{  4  D  )  / 


(291 


Table  1  compares  values  of  2 Q/kxa  for  the  various  predictions  (27)-(29)  while 
Fig.  1  plots  the  percentage  deviation  between  the  predictions  (29)  and  (27).  It  can  be 


o  J/Dt  or  a  JWj/0 

Fig.  I.  Hoc  of  the  retail**  error*  with  re* pact  to  th«  eua  CE  result  (27)  or  the  exact  cfaropomperometnc 
result  (?)  (or  (i)  tbs  CE  approximation  in  eqn.  (29)  and  (b)  (or  the  ebrooonmpcrometjic  approximation 
in  eqn.  (33). 


seen  that  the  approximate  method  gives  results  which  are  in  fact  very  close  to  the 
result  (21). 

CONCLUSION 

The  present  series  of  papers  (see  also  refs.  3-5)  has  shown  that  it  is  possible  to 
analyze  a  wide  range  of  electrochemical  experiments  at  microdisk  electrodes  by 
making  use  of  the  properties  of  discontinuous  integrals  of  Bessel  functions.  In  this 
analysis  we  have  made  extensive  use  of  the  average  concentration  over  the  surface  of 
the  microdisk  rather  than  solving  the  problems  for  controlled  uniform  concentra¬ 
tions  over  the  surface.  The  latter  approach  (although  possible)  leads  to  considerable 
mathematical  complications  and  results  which  frequently  have  to  be  approximated 
It  is  of  interest  ihuitfore,  that  the  approximate  procedure  used  previously  for  the 
discussion  of  coupled  reactions  in  solution  (in  which  the  exact  result  for  me  problem 
for  a  microspherical  electrode  is  scaled  to  that  for  the  disk  using  eqn.  (3)  [3])  gnes 
results  which  are  close  to  those  using  the  procedure  developed  in  this  senes  of 
papers.  As  a  further  illustration  of  this  procedure,  we  note  that  the  diffusion  limited 
current  density  to  a  sphere  (a  special  case  of  constant  surface  concentration) 

I  -  nFDc*/a^„t  +  nFDu/ic* /nl/2tl/2  <  " 

can  be  rescaled  to  that  at  the  disk  using  eqn.  (3). 

/  -  4 nFDc* /irailik  +  nFD'/2c~/ir'/2tl/ 1  l  ' . 


8 

We  therefore  obtain  the  actual  current  at  a  disk  of  radius  adlik  as 

ira 3„k/  -  4 nFDc^a^  -r  nFir1  :D1'  2 

The  current  transient  normalized  by  the  steady  state  current  is 

_  zLL  i  j 

4  nFDe*  “  "  4  'o'1'2!*'2' 


(32 1 


133) 


We  note  that  eqn.  (33)  is  close  to  the  values  predicted  algebraically  [6]  or  bv 
simulation  (7]  the  maximum  deviation  being  +6.8%  at  aiak/Dl/1tx‘ 2  a  3  this 
devtauon  decreasing  at  both  high  and  low  values  of  a/D'/1t'/1.  Table  1  and  Fig.  1. 
We  note  that  a  wide  range  of  problems  has  already  been  analyzed  in  the  spherical 
coordinate  system  and  that,  when  solutions  are  not  available,  these  are  easy  to 
obtain  in  view  of  the  uniform  accessibility  of  the  surface  of  spherical  electrodes.  The 
approximate  procedure  may  therefore  prove  to  be  generally  useful. 
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